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We utilize caustic theory in "PT—symmetric lattices to design focusing and curved beam dynamics. 
We show that the gain and loss parameter in these systems provides an addition degree of freedom 
which allows for the design of the same caustics trajectories with different intensity distribution 
in the individual waveguides. Moreover we can create aberration-free focal points at any paraxial 
distance Zf, with anomalously large focal intensity. 


Introduction- Diffraction management of beam propa¬ 
gation and the possibility of designing initial wavefronts 
which lead to abrupt energy focusing have attracted the 
attention of many researchers over the years [I]j8]. Not 
only is the fundamental and mathematical physics side 
of this problem both charming and challenging for re¬ 
searchers, but also its applied side has attracted atten¬ 
tion. Perhaps the most pronounced example comes from 
the field of medical lasers where one needs to have an 
abrupt beam focusing at a specific point, without affect¬ 
ing nearby tissues [3]. Other applications include parti¬ 
cle manipulation |10j . generation of self-bending plasma 
channels m and optical bullets [HIS] etc. 


n-1 n n+i 



K C K C K 


FIG. 1. A PT-symmetric lattice consisting of dimers with 
local PT-symmetry. Each dimer consists of two waveguides, 
one attenuating (type ”a” waveguides, indicated with green) 
while the other one (type ”b” waveguides, indicated with 
red) has equivalent amplification. The coupling between two 
waveguides of a specific dimer is k while the coupling between 
different dimers is c < k. 

In this Letter we will investigate caustics trajectory 
dynamics and abrupt focusing in parity-time (VT)- 
symmetric discrete array settings, like the one shown in 
Fig.[T} These type of optical systems have been recently 
introduced [Hj and, during the last years, have gained a 
lot of attention due to the wealth of exotic properties that 
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they posses [IlHSllEelETlESHaS]. The main characteris¬ 
tic of these systems is that the optical potential (index of 
refraction) e{x) = eji{x) ij{x) is complex and satisfies 
the PT symmetric property e{x) = e*{—x). Since the 
optical potential is complex these systems are described 
by an effective non-Hermitian Hamiltonian with eigen¬ 
values (propagation constants in optics language) which 
are real if the gain/loss parameter is below some critical 
value i.e. 7 < uvT- Above this value the eigenvalues be¬ 
come complex and the system is unstable. The transition 
point 7 = j-pT has the properties of an exceptional point 
(EP) singularity i.e. both eigenvalues and eigenvectors 
coalesce. This phase transition (coined PT -symmetric 
phase transition) leads to a number of interesting fea¬ 
tures: asymmetric transmission [nHu, unidirectional 
invisibility [22U24] , novel lasing schemes [25l [26l [28] , non¬ 
reciprocal Bloch oscillations [331 El] j and reconfigurable 
Talbot effects [35] . 

In this paper we will demonstrate the possibility of de¬ 
signing caustics and identify the effects of non-hermiticity 
(i.e. gain/loss elements described by a complex index of 
refraction) in engineering beam trajectories and in the 
formation of abrupt focusing processes. Our approach 
will utilize curved trajectory dynamics and caustics de¬ 
sign in discrete elements and thus extend the analysis for 
passive lattices [5] to lattices with PT -symmetry. 

Model and theoretical analysis- We consider an array 
of one-dimensional (ID) coupled waveguides. The array 
is composed of two types of waveguides: the first with 
gain (A) and the other with equivalent loss (B). We fur¬ 
ther assume that each waveguide is supporting only one 
propagating mode. The nearby waveguides are assumed 
to be coupled evanescently. The (A) and (B) waveguides 
are arranged in a way that they form N coupled A-B 
dimers with intra-dimer coupling k and inter-dimer cou¬ 
pling c. In the slowly varying envelope approximation the 
electric field amplitude = (a„, bn)'^ at the n-th dimer 
evolves (along the propagation direction z) according to 
the Schrodinger-like equation: 

(z) kbn{z) + cb„_i(z) = 0 

^ db„T) _|_ Cbn{z) + kan{z) + can+l{z) = 0 

where the complex refractive index of the n-th waveguide 
is e = eo + *7. eo is the background index of refraction 
and 7 is the gain-loss parameter. Without any loss of 
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generality, we will assume below that eo = 0. Equation 
Q , which describes the field evolution, is invariant under 
a parity-time {VT) symmetric operation [3T1 [35] . The 
parity-symmetry operator V is defined as a spatial inver¬ 
sion (around an axis of symmetry of the array) while the 
time-reversal symmetry operator 7" is associated with a 
complex conjugation. 

The propagation of the electric field along the parax¬ 
ial direction z is conveniently evaluated in the Fourier q 
space where a„(z) = ^ dqdg(z) exp(ing) (similarly 
for bn). Specifically, the translational invariance of the 
system allows us to decouple, for each value of q, the 
equations of motion into 2 x 2 blocks: 


•7^ / dq{z) 
^dz V hq{z) 



\H„ 


-*7 Vq \ 

V* Z7 J 


( 2 ) 


where Vq = —(A:-|-c-e“®*). Equation ([^ can be solved an¬ 
alytically, thus allowing us to evaluate the wave packet 
evolution in the Fourier space. The wavepacket in the 
real space is then obtained by an inverse Fourier trans¬ 
form dq{z) = J2m=-oo (similarly for bq{z)). 

The dispersion relation of the dimeric lattice is calcu¬ 
lated from Eq. ([^ by substituting the stationary form 
{an,bn)^ = exp{—i£z){A,B)'^ [ST]. We get: 


^±( 9 ) = i \/{k — c)^ -I- 4:kccos'^{q/2) — 7 ^. (3) 

It follows from Eq. ([^ that for 7 < 7737 - = k — c the 
dispersion relation consists of two bands which are sepa¬ 
rated by a gap. In this parameter domain all the eigen¬ 
values are real and the system is stable. The maximum 
gap size 2{k — c) occurs for 7 = 0. For larger values of 7 
the gap becomes smaller until it disappears at ^vT- At 
7 = IVT fhe levels associated with q = ± 7 r and their 
corresponding eigenvectors become degenerate resulting 
in an exceptional point (EP) singularity. For 7 > j-pT 
the spectrum becomes partially complex [81j . Bellow we 
assume that k > c. 

Using the dispersion relation Eq.(|^ we evaluate the 
electric field amplitudes associated with the n—th dimer 
'kn = {an(z),bn{z))’^ at any distance z: 


\buiz)) 



iq{n—m) — iS±z 


dq. 


(4) 

To proceed with our analysis we extend the integer index 
TO, defining the dimer number, to a continuous variable 
^ which can be interpreted as the transverse spatial co¬ 
ordinate variable of the incident (i.e. z = 0) field. Simi¬ 
larly, we express the integer index n as a general continu¬ 
ous transverse spatial-coordinate x. This extension allow 
us to define the continuous smooth functions a(^, 0 ) and 
6 (^, 0 ) with the property that a(|^| —>■ oo, 0 ) —>■ 0 and 
similarly 6 (|^| —>• oo,0) —0. Next we express Eq.Q in 
the following integral form: 


b{x,z))-2nJ_^J_A^d 

( 5 ) 


where we have used the polar representation a(^, 0 ) = 
and 6(^,0) = The subindex ^ means that 

a, j3 and (j) are functions of It turns out that the caus¬ 
tics formation is independent of the amplitudes a^; (3^ 
(see Eqs. below) and therefore we will assume 

below that aj = = 1 . 

In order to enforce caustic beam dynamics, we impose 
a stationary phase (first and second order) condition on 
both integration variables q and ^ in Eq. ([^: 


d<s> 

di 


^=0^q^ = ^ (a) 

H = 0 ^a; = e-g^sinq^ {b) 

Along with the second order stationary condition 


( 6 ) 




_= 0 ^ Z = _ 

dq'^ ^ dq ’ 2(c^-t-fc^—7^) cos q^+cfc(3-t-cos 2g^) 

(7) 

where q'^ = ^ and $ = -I- q{x — — £±z. Equations 

define the coordinates x and z of the ray trajectory 
associated with the ^th dimer. They allow us to design 
various types of caustics associated with the Eqs. Q of 
a beam propagating through the UT-symmetric dimeric 
lattice of Fig. Solving Eq. ([§3) with respect to q^ 
leads to the expression: 


where rj = . Further integration of Eq. (j^ with re¬ 

spect to ^ (see Eq. ([^)) results in a relationship between 
the lattice parameters, the ray trajectories, and the ini¬ 
tial phases of the waveform cj)^. This enables us to tailor 
the initial phases cj)^ in order to achieve a desired caustics 
dynamics. The condition that the initial phases must 
be real imposes a constraint to q^ that allows for the cal¬ 
culation of the maximum number of dimers ^ < |^o| that 
participate in the caustic dynamics. 

We start with designing a focusing point at a prede¬ 
fined transverse coordinate position x / and at a paraxial 
distance z/. Below we keep z/ as a free parameter while, 
for the sake of presentation, we assert that Xf = 0, i.e. 
it coincides with the origin of the transverse coordinate. 
The requirement for focusing at a;/ =0 further simplifies 
Eq. (|^, which can be used to evaluate the initial phase 
of the wavefronts. We get: 

4=±|i|£±(97)|±Cg/ + z/Co (9) 


f — 1 / 2 

= cos 


[ (c 2 -h - 72),;2 ^ 

V -s—‘ + d) 


( 10 ) 


where Co = a/(c -I- k)^ — 7 ^ and r]f = . The ex¬ 

istence of a focal point for z/ > 0 is achieved when the 
initial phases are for ^ < 0(^ > 0). The number 
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of dimers |^o| that participate in the creation of the focal 
point can be calculated from Eq. Q. We get 


+ fc2 - 72 - Co v^r - 7 ^ 

Co - ---- (11) 



FIG. 2. Design of an aberration-free focal point at Zf = 20 
of an initial wavefront for a 'PT-symmetric lattice with two 
different 7-values. In (a) and (c) the intensity profile of the 
evolving engineered beam is presented. In (b) and (d) we 
show the associated ray trajectories. Each individual ray cor¬ 
responds to a dimer. The initial wavefront in (a) and (c) has 
the form a(^,0) = and b(^,0) = for Co < C < |Co|, 
where the initial phases are given from Eq.H. The cou¬ 
pling constants are k = 5, and c = 1. In (a) and (b), 7 = 0 and 
therefore we use |Co| ~ 20. In (c) and (d), 7 = 3.999 ~ 'yvT 
and hence |Co| = 44. In this case the field intensity at the 
focal point is an order of magnitude larger than in (a) 


Equations ( |9|10|11| allow us to gain a better under¬ 
standing of the role of the gain and loss parameter 7 in 
the formation of focusing points. Let us consider, for 
example, the two scenarios of a case with 7 = 0 and a 
VT- symmetric case with 7 = 'y-pT- In fh® former case 
we find from Eq. (II) that |Co| = czf dimer trajectories 
(and thus dimer waveguides) contribute for the construc¬ 
tion of a focal point at a paraxial distance Zf from the 
plane of preparation. In contrast when 7 —> 'j-pj- we 
gst ICol = 's/cfcz/ i.e. more waveguides are contributing 
to a focal point. In Fig. we show the evolution of 
an initial beam which is designed to create an ideal (i.e. 
aberration free) focal point for the two different values 
of 7 . The phase engineered of the initial beam is dic¬ 
tated by the rules given by Eqs. ([^ |I0|II ) for each 7 
value. We find that increasing the gain/loss parameter 
7 leads to an increase of the number of waveguides ^0 
which contribute to the creation of the focal point. At 
the same time the total beam intensity at the focal point 
increases dramatically. This growth is a consequence of 
the EP dynamics and has been predicted theoretically in 
m and observed experimentally in |24j . In the current 
framework this growth is manifested as a giant intensity 


growth at a pre-engineered focal point. 



FIG. 3. Two representative examples of designed caustic 
propagation juxtaposed with the associated caustic trajectory 
of the form x = az^ . In both cases the lattice parameters are 
kept fixed fc = 3, c = 1, and 7 = 1.95. The initial condition 
for the system is a(^, 0) = b(^, 0) = exp"^? for ^0 < ^ < 0 
and zero elsewhere. The difference between (a) and (c) is the 
caustic trajectory used to calculate the initial phase 0^. (b) 
is the associated ray trajectory for (a) and (d) is the asso¬ 
ciated ray trajectory for (c). In (a) d = 3, a = 0.0003 and 
consequently ^0 ~ —40, while in (c) <5 = 5/4, a = 0.4 and 
consequently ^0 ~ —25. The green lines in (b) and (d) are 
the desired caustic trajectories. 


Next we investigate the management of curved caus¬ 
tic trajectories in T^T-symmetric lattices. We start with 
Eqs. ( |6|7|8 ). For theoretical simplicity the ensuing calcu¬ 
lation will focus on designing caustics associated with an 
initial wave-front involving only dimers with index ^ < 0 . 
We will assume a power law caustic trajectory x = az^. 
Furthermore we relate the caustic propagation distance 


z, with the dimer index ^ using Eq. &) which can be 
re-written in a more compact form as ^ = x — z^. One 
then uses these two relations to write the variables x and 
2 : as functions of the dimer index the power S that de¬ 
fines the caustic trajectory, and the scaling constant a. 
The resulting expressions are: 


X = 


(I-h)’ 


= (- 


:(I-5) 


)l/« 


( 12 ) 


The above relations encode the information of the caus¬ 
tic trajectory into the ray trajectory coordinates for the 
^-th dimer. Substituting these expressions back into Eq. 
([^, enables us to calculate the initial phase distribution 
producing dynamics following the caustic trajectory 
defined by Eqs. (12). Furthermore using Eq. we 
estimate the number of dimers ^0 that are used to con¬ 
struct the propagating caustics. In Fig. we present two 
examples of different caustic trajectories for two differ¬ 
ent S values. Our results confirm that power-law caustic 
trajectories can be supported in T^T-symmetric lattices. 
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Obviously this agreement is valid up to some paraxial 
propagation distance z and transverse coordinate vari¬ 
able X which are given by Eq. ([T^ for ^ 


We can also study the effect of the gain/loss parame¬ 
ter 7 on a specific caustic trajectory. This is presented 
in Fig 1^ where we demonstrate the beam dynamics orig¬ 
inated by the same caustic trajectory for two different 
values of gain/loss parameter. We see that increase of 
7 , while also adjusting the initial phase accordingly, aug¬ 
ments the amplitude of the caustic wave without altering 
its trajectory. 


FIG. 4. Two examples of designed caustic propagation juxta¬ 
posed with the desired caustic trajectory of the form x = az^, 
where a = 0.1 and & = 3/2, for (a) 7 = 0.5 and (c) 

7 = 1.8. In both cases the lattice coupling parameters are 

kept fixed and equal to fe = 3, c = 1. The initial wavefront is 
0 ) — 0 ) = exp*"^? for ^0 < C < 0 a^nd zero elsewhere. 

In (a) ^0 = —15 while in (c) ^0 = —36. (b) is the associated 

ray trajectory for (a) while (d) is the associated trajectory 
for (c). In (a)-(b) zjg = 45 and = 30 while in (c)-(d) 
zjo = 80 and x^g = 72. Furthermore, the point of maximal 
intensity along the caustic in (c) occurs approximately at the 
same coordinates as the end of the caustic in (a). The green 
lines in (b,d) are the desired caustic trajectories. 


In conclusion we have investigated the possibility of 
implementing caustics design in T^T-symmetric lattices. 
We have demonstrated the design of power law caustic 
trajectories and the possibility of creating aberration-free 
focal points with increased focal power due to the pres¬ 
ence of EP singularities. 
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